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Abstract 

A program including all radiative corrections to the Minimal Supersymmetric Standard 
Model (MSSM) at the same level as the radiative corrections to the Standard Model (SM) 
has been developed and used to perform global fits to all electroweak data from LEP, SLC 
and the Tevatron and the radiative b — > 57 decay from CLEO. Values of the strong coupling 
constant at the Mz scale and sin 2 9jjg are derived, both in the SM and MSSM. 

Recent updates on electroweak data, which have been presented at the Warsaw Confer- 
ence in summer 1996, reduce the too high branching ratio of the Z° boson into b-quarks 
in the SM from a 3.2a to a 1.8a effect. In addition, the b — > s'j decay is 30% below the 
SM prediction. In the MSSM light stops and light charginos increase Rb and decrease the 
b — > 57 rate, so both observations can be brought into agreement with the MSSM for the 
same region of parameter space. However, the resulting y 2 value for the MSSM fits is only 
marginally lower and in addition, the splitting in the stop sector has to be unnaturally high. 
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1 Introduction 



Super symmetry presupposes a symmetry between fermions and bosons, which can be im- 
plemented in the Standard Model (SM) by introducing a fermion for each boson and vice 
versa 0. In this case the problem of quadratic divergent radiative corrections to the Higgs 
boson masses is solved, since fermions and bosons contribute with an opposite sign to the 
loop corrections. These new supersymmetric particles ( "sparticles" ) contribute addition- 
ally to the radiative corrections and can influence electroweak precison variables, like e.g. 
Rb = r z°^bi/^ ' z°->hadrons or the penguin mediated decay b — > sj. Radiative corrections have 
been calculated in the Minimal Supersymmetric Standard Model (MSSM) to nearly the same 
level as in the SM, so an equivalent analysis of all electroweak data can be performed both 
in the SM and MSSM. In this paper such analysis are described using data from Tevatron 
@, |], LEP and SLC @, |, the measurement of R b ^ sl = f^^| from CLEO § and 
limits on the masses of supersymmetric particles [7L || |J ITDL 11]. For the SM predictions 



the ZFITTER program[12fl was used, while for the MSSM predictions the SUSYFITTER 



program||13||, which will be discussed below, was used. In both cases MINUITJH1 was used 



as x minimizer in order to obtain the optimum parameter values. For all SUSY masses well 
above the electroweak scale one does not expect significant differences between the SM and 
the MSSM predictions. 

Previous LEP data || showed a too high value of Rb (3.2a) and a too low value of R c 
(2a). It has been shown by several groups [[OJ-|23) that it is possible to increase Rb in the 



MSSM with light charginos, top squarks or Higgses, which yield large positive contributions 
to the Zbb vertex because of the large Yukawa couplings of the third generation. The 
first two generations are not affected by such corrections, so no modifications in R c can be 
obtained. Recent updates of electroweak data, presented at the Warsaw Conference B], show 
no deviation of R c anymore and a value of Rb which is 1.8a above the SM value. In addition, 
the experimental value for i2&_» S7 = §§^E^k from CLEO U equals (2.32±0.67) ■ 1CT 4 , which 
is about 30% below the SM prediction after taking the ca. 10% increase in the prediction 



by the next-to-leading-log contributions into account 1 24}] . 

In the MSSM an increase in Rb can cause a decrease in the b — > 37 rate, since both cases 
involve similar diagrams usually with an opposite sign, e.g. the t — t — x =l= vertex corrections 
in Z° — > bb and the t — X s1 loop corrections in b — > 57 have an opposite sign. Both, the too 
high value of Rb and the too low value of b — > S7 prefer a chargino mass around 90 GeV and a 
stop mass around 50 GeV or alternatively light Higgses around 50 GeV, so both observations 
agree with the MSSM for the same region of parameter space, albeit an unnatural one, as 
will be shown. 

2 Radiative Corrections in the MSSM 

At the Z boson resonance two classes of precision observables are available: 
a) inclusive quantities: 

• the partial leptonic and hadronic decay width T^j, 

• the total decay width T z , 
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• the hadronic peak cross section 0^, 

• the ratio of the hadronic to the electronic decay width of the Z boson: R^, 

• the ratio of the partial decay width for Z — > cc(bb) to the hadronic width, R c , 

Rb- 

b) asymmetries and the corresponding mixing angles: 

• the forward-backward asymmetries A FB , 

• the left-right asymmetries A LR , 

• the t polarization P T , 

This set of precision observables together with the Rb^ s -y rate|J, the W - and the top mass is 
convenient for a numerical analysis of the supersymmetric parameter space. In the following 
the observables defined above are expressed with the help of effective couplings. 



2.1 The effective Z-f-f couplings 



The coupling of the Z boson to fermions / can be expressed by effective vector and axial 
vector coupling constants ^e//> a eff m terms of the NC vertex: 



(1) 

I M\ |25|. Input 

parameters are the \i decay constant = 1.166392 x 10~ 5 GeV -2 , ocem — 1/137.036 and 
the mass of the Z° boson Mz = 91.1884 GeV. The mass of the W boson is related to these 
input parameters through: 



where the convention is introduced : c™ 
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cos 2 6 



' w 



Gjj_ _ KOtEM 1 

y/2 2s^M^ 1 - Ar MSS M (a em , M w ,M Zl m t ,...) 



(2) 



where the complete MSSM one-loop contributions are parameterized by the quantity 
ArMSSM • Leading higher order SM corrections |27|, |28| to the quantity Ar are included 
in the calculation. 

The effective couplings v{ 



j e jj-, a{,ff can be written as: 



J eff 



'Z z (v f + Av f + Z M Q f ) 

"eff - V Zz + ' 

v* and a* are the tree-level vector and axial vector couplings: 



I( - 2Q f s 2 w ,af = l 



f 
3 • 



(3) 



(4) 



Z z , Z M are given in eq. (fL0|). The complete MSSM one- loop contributions of the non- 
universal finite vector and axial vector couplings Av? , Aa? have been calculated JT3|, together 
with the leading two-loop Standard Model contributions [27, 28, 2^]. They are derived in the 



't Hooft-Feynman gauge and in the on-shell renormalization scheme [pO 
MSSM one-loop Z — > ff vertex correction diagrams. 



Fig. 1 shows the 
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Figure 1: MSSM one-loop Z — > // vertex correction diagrams. i,j, k = 1, ..,2(4) are 
chargino (x*), neutralino (x°) and sfermions (/) indices, while h, H, are the scalar Higgs 
bosons, A is the pseudoscalar Higgs boson, and G are Goldstone bosons. No particle per- 
mutations are shown. , 



y,Z° W + h°,H ,A ,G' 

rJ^M^ rJ^^U --"-^ 



fff fff fff 




The non-universal contributions can be written in the following way: 



Av f = Fy M + AFy, 
Aa f = F% M + AF A . 

The Standard Model form factors Fy 1 ^ corresponding to the diagrams of figs. 1 and 2 can be 
found e.g in refs. [[H], |2"7| . The diagrams with a virtual photon are listed for completeness in 
the figures. They are not part of the effective weak couplings but are treated separately in the 
QED corrections, together with real photon bremsstrahlung. The non-standard contributions 
are summarized by 

AFy = Yl F v +v f SZ v + a f 6Z f A 

i 

AF A = J2 F a +v f 5Z{ + a f 5Z v 

i 

where the sum extends over the diagrams of fig. 1 with internal charged and neutral Higgs 
bosons, charginos, neutralinos and scalar fermions. Each diagram contributes 

Fv% - FA 7m75 

to the Z f f '-vertex. The self-energy diagrams of fig. 2 with internal neutral Higgs, chargino, 
neutralino and sfermion lines determine the field renormalization constants 



5Z f A = £ A K) (5) 



5Z V = -E v (m 2 )-2m}[E' v (m}) + E' s (m})} 
with the scalar functions £\r,A,s m the decomposition of the fermion self-energy according to 

The contributions from the Higgs sector are given explicitly in ref. f3]]. For the genuine 



SUSY diagrams, the couplings for charginos, neutralinos and sfermions are taken from [[32 
together with the diagonalization matrices given in the next section. 



The universal propagator corrections from the finite Z boson wave function renormaliza- 
tion Zz and the 7Z mixing Zm are derived from the (7, Z) propagator matrix. The inverse 
matrix is: 

' k 2 + t 1 {k 2 ) t lZ (k 2 ' 
±,z(k 2 ) k 2 -M 2 +^ z (k 2 



(vr i =<^ izzr* ,r\L , *m , (6) 



where S 7 , S^, are the renormalized self energies and mixing. They are obtained by 
summing the loop diagrams, shown symbolically in fig. |3|, and the counter terms and can be 
found in ref. [I3|. 



The entries in the (7, Z) propagator matrix: 
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are given by: 



A z (k 2 ) 
A, z (k 2 ) 



k 2 + S 7 (P) - 



fc 2 -M 2 +E z (fc 2 ) 



k 2 - M 2 Z + ± z (k 2 



fc 2 +E 7 (fc 2 ) 

£ lZ (k 2 ) 



[k 2 + E 7 (k 2 )} [k 2 - Ml + E z (k 2 )} - £ 2 z (fc 2 ) 



(8) 



The renormalization condition to define the mass of the Z boson is given by the pole of the 
propagator matrix (eq. []). The pole k 2 = Mf is the solution of the equation: 

Tie [ (M 2 Z + S 7 (Ml) ) £ z (k 2 ) - %(Ml) } = . (9) 

Eq. (R|) yields the wave function renormalization Z z and mixing Zm'- 



Res Mz A z 



1 



1 + Tle£' z (k 2 ) - lie 



^z(k 2 ) 



k 2 +E 7 (k 2 ) 



fc 2 =M§ 



z 



%z{M 2 ) 



M 2 + S 7 (M|) 
2.2 Z boson observables 

The fermionic Z boson partial decay widths V can be written as follows: 
1) / ^ b: 



(10) 



where 



2) / = b: 



ff 



NcVlG^M* 



127T 

3o>em 
An 



1 - Atmssm) 



v f eff ) 2 + (al„) 2 (l- 



effJ 



6mj 
Ml 



Q 2 f ) (i + *£ 



qcd; 5 



S f 

°QCD 



, / = leptons 

^ + 1.405(^) 2 - 12.8(^) 3 - , / = quarks 



NcV^G^Ml 



12tt 



(1 - Ar MSSM ) 



+ (a b eff ) 2 



{l + ^Ql){l + 8 b QCD )+AT bl . 



(11) 
(12) 



(13) 



In AT b i the b quark specific finite mass terms with QCD corrections are included. 
Sqcd is g iven in e Q- ®- 



The total decay width T z is the sum of the contributions from leptons and quarks: 

r ^ = E r //- (14) 

/ 

In the following Thad = J2 q F qq - is the hadronic decay width of the Z boson. 
The hadronic peak cross section is defined as 

Yin T ee Y had 

° h = Mi-rr ■ (15) 

The ratio of the hadronic to the electronic decay width is defined as 

Re = ^ • (16) 
J- ee 

The ratio of the partial decay width for Z — > bb (cc) to the total hadronic decay width is 
given by 

RHc) = ^ • (17) 

J- had 

The following quantities and observables depend on the ratio of the vector to axial vector 
coupling. The effective flavour dependent weak mixing angle can be written as 

/ 



a eff, 



4IQ/I 

The left-right asymmetries are given by 

jf eff/ a i 

1 + (vlf f /al ff y ' 

while the forward-backward asymmetries can be written as 



Af - /If- 2v eff/ af eff 

^-lr — • /i _ " Ti ; f 



(19) 



4 B = ^M f . (20) 



3 The MSSM 

3.1 Higgs sector 

The MSSM has two Higgs doublets: 



_ ( v i + 75 (H° cos a - h° sin a + iA° sin (3 - iG° sin (3) \ 
1 ~\ ' H~ sin (3 - G~ cos (3 J ^ ' 

( tf+cos^ + G+sin/3 \ 

n2 ~ y V2 + -± (H° sin a + h° cos a + iA° cos (3 + iG° sin P) J {ZZ) 



Here H, h and A represent the neutral Higgs bosons, while the G"s represent the Goldstone 
fields, which correspond to the longitudinal polarization components of the heavy gauge 
bosons. The imaginary and real sectors do not mix, since they have different CP-eigenvalues; 
a and (3 are the mixing angles in these different sectors, so one is left with 2 neutral CP- 
even Higgs bosons H° and h°, one CP-odd neutral Higgs bosons A , and two charged Higgs 
bosons. 

Their masses are completely determined by the ratio of the vacuum expectation values of 
tan/3 = V2/V1 and the pseudoscalar mass Ma, together with the radiative corrections. The 
latter ones are taken into account in terms of the effective potential approximation with the 
leading terms ~ mf, including the mixing in the scalar top system |j33"| . In this way, the 



coupling constants of the various Higgs particles to gauge bosons and fermions can be taken 
over from [52] substituting only the scalar mixing angle a by the improved effective mixing 
angle which is obtained from the diagonalization of the scalar mass matrix, discussed in the 
next section. 



3.2 Sfermion sector 

The physical masses of squarks and sleptons are given by the eigenvalues of the 2x2 mass 
matrix: 

V| , , M 2 Q + m} + MUll-Q f s 2 w )cos2f3 m f (A f + p{cot 0, tan/3}) 
" > 1 mf(A f + fi,{cotp,tmP}) M 2 {(/ D} + m) + M 2 z Q f s 2 w cos 2/5 » ' ' 



with SUSY soft breaking parameters Mq, Mjy, Mp, Af, and the mass parameter \i from the 
Higgs sector It is convenient to use the following notation for the off-diagonal entries in 
eq- ©: 

A' f = A f + p{cot P, tan 0} . (24) 

Scalar neutrinos appear only as left-handed mass eigenstates. Up and down type sfermions 
in (p3|) are distinguished by setting f = u (d) and select cot /3 (tan/3) in the curly brackets. 
Since the non-diagonal terms are proportional to rrif, it seems natural to assume unmixed 
sfermions for the lepton and quark case except for the scalar top sector. The t mass matrix 
is diagonalized by a rotation matrix with a mixing angle $> m i X . Instead of Mq, My, M5, A' t 
for the b, t system the physical squark masses m^,m^, m f ~ 2 can be used together with A' t 
or, alternatively, the stop mixing angle $ m i X . For simplicity we assume m~ bL = rn~ bR , and u, 
d, c, s to have masses equal to the b squark mass. 

A possible mass splitting between b^-t^ yields a contribution to the p-parameteiQ p = 
1 + Ap° in terms of f||: 



a 3a EM (2,2 o b L m t L , m b L \ rn ,s 

r = tt. — 9 ,„ rn~ + m~ - 2^— L - — V log . (25) 

6 ~* l&ns 2 w M^ \ bL tL m| - m\ h mj L J v ' 

As a universal loop contribution, it enters the quantity 

Ar ~ Aa EM -^fAp° + ... (26) 
s w 



°The superscript in Ap° indicates that no left-right mixing is present. 
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and all the Z boson widths 



r ff -~ i + Ap° + 



and is thus significantly constrained by the data on M w and the leptonic widths. 



3.3 Chargino/Neutralino sector 

The chargino (neutralino) masses and the mixing angles in the gaugino couplings are cal- 
culated from jj, and the soft breaking parameters Mi, M 2 in the chargino (neutralino) mass 
matrix p2]. 

The chargino 2x2 mass matrix is given by 



M 2 M w V2sm(3 
MwV^ cos (3 —fi 



(27) 



The physical chargino mass states xt are the rotated wino and charged Higgsino states: 



xt = Viiift 

Xi = UijipJ ; i, j = 1,2 . 



(2? 



Vij and Uij are unitary chargino mixing matrices obtained from the diagonalization of the 
mass matrix eq. ^7|: 

ITAVV- 1 = diag(m-±, m-±) . (29) 



The neutralino 4x4 mass matrix can be written as: 



Mi -M z sin W cos (3 M z smO w sm(3 \ 

M 2 M z cos 9w cos (3 —M z cos 6 W sin (3 

-Mz sin 6w cos (3 Mz cos Ow cos (3 \i 

\ Mz sin Ow sin (3 —Mz cos Ow sin (3 \i / 

(30) 

where the diagonalization can be obtained by the unitary matrix Nif. 

N^^oN" 1 = diag(m^o) . (31) 

The elements U^, Vij, Ny of the diagonalization matrices enter the couplings of the charginos, 



neutralinos and sfermions to fermions and gauge bosons, as explicitly given in ref. [32|. Note 



that our sign convention on the parameter /i is opposite to that of ref. [[32 . 



4 Results 



The experimental limits included in the fit are summarized in table |l] Jf[ §, || |10|, |TT|. They 
have not been updated with the latest results, since in practice only i4->s7 an d Rb require 
light SUSY particles, but their deviation from the SM has become so small, that they do 
not require sparticles below the experimental limits. The experimental observables used in 
the fits are shown in the first column of table 0. The calculation of the total decay width of 
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experimental limits 


m-± 

Al,2 


> 65 GeV 


A3, 4 

Tz— >neutralinos 


> 13 GeV 

> 35 GeV 

> 60 GeV 
< 2 MeV 




> 48 GeV 


m h ,m H ,m A ,m H ± 


> 50 GeV 



Table 1: Mass limits assumed for the optimized fits. 



34 , the calculation of the ratio Rh 



the Z boson into neutralinos is based on reference 
reference j35|. The next-to- leading- log calculations increase the SM prediction for i?6_» S7 by 
about 10% to (3.2±0.5)10~ 4 |p4| . This higher order contribution was taken into account in 
the first order calculation by choosing the renormalization scale fx = 0.7 • m&. 

As discussed in section p4j , Rb depends on stop -, chargino- and Higgs masses. First 
the behaviour of the chargino masses as a function of the SUSY parameters fx and M 2 is 
discussed, then the Rb dependence on the relevant SUSY masses is studied, both for the low 
and high tan/3 scenarios. These studies are followed by the best solutions in the SM and 
MSSM. 



4.1 Chargino Masses 

In order to explain an enhanced value of Rb in the data, the MSSM needs a light right 
handed stop and light chargino (low tan/3 scenario), or a light pseudoscalar Higgs A (high 
tan/3 scenario) |]15||-||23||. A higgsino-like chargino can be obtained for a low value of the 



parameter fx in the mass matrix (eq. [27]). Figs. [| and |5] show the dependence of the 
chargino masses on the parameter fi. For high tan/?, m^ 2 is almost symmetric around fi = 0, 
whereas for low tan /3 this dependence is more complicated, as can be seen from fig. f|. For 
Mi = 3\ix\ the light chargino mass passes zero at fx — —40 GeV, so the following low tan/3 
plots were made for fi > —40 GeV and fi < —40 GeV. The asymmetric structure of fig. f| 
is reflected in the contours of constant Rb in the m^ 2 versus light scalar top m f ~ 2 plane (see 
fig. |S|). Values of Rb up to 0.2194 are possible (see figs. |G] and [7p, although these special 
regions of the parameter space are already experimentally excluded by the lower limits on 
sparticle masses. 

For M 2 = 3\fx\ the lightest chargino is mostly higgsino-like, while the heavier one is 
gaugino like. Mixing the charginos more by taking M 2 = \fi\ does not change these results 
very much, as can be seen from a comparison of the x 2 distributions in fig. |8| (M 2 = 3|/x|) 
and fig. P (M 2 = \fi\). The small increase of the x 2 a ^ chargino masses around 80 GeV in the 
left hand part of fig. |9] is due to neutralino threshold singularities, for which an additional 
X 2 contribution has been added, when the sum of two neutralino masses is close to the Z° 
mass. The sharp increase of the x 2 function at low chargino masses is due to experimental 
limits on chargino, neutralino and stop masses from LEP 1.5 0, §, [J. 
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4.2 Low tan/3 scenario 

Fig. [10] shows the change in the best obtainable \ 2 in t ne chargino - stop plane. For each 
value of the lighter scalar top m t - 2 and lighter chargino m~± in a grid of 10x10 points an 
optimization of m t , a s and the stop mixing angle $ m j X was performed, assuming M 2 = 3\/i\ 
for a fixed value of tan/3 =1.6. In the next section this assumption on M 2 will be dropped. 
Low sparticle masses yield a sharp increase in the A% 2 in fig. pi] because of the included mass 
limits. The minimum x 2 is obtained for chargino masses above 80 GeV and increases only 
slowly with increasing sparticle masses. Rb increases significantly with decreasing values of 
the stop and chargino mass, as can be seen from fig. [11]. Much less significant is the change 
of R c . Within the plane of fig. O it changes less than 0.0005 units. The increase of Rb 



must be compensated by a decrease of a s (see fig. O) in order to keep the total Z°-width 



constant. The stop mixing angle & m i X , shown in fig. [13], is mainly determined by the CLEO 
measurement of Rb-^sj- The chargino contribution to Rb-> s ~/ is proportional to the Higgs 
mixing parameter /i, which changes its sign for m^± ps 60 GeV (see fig.fD, so the Rb^ s ^ rate 
changes rapidly for these chargino masses, as shown in fig. |14] 



4.3 High tan/3 scenario 

Similar fits can be performed for the high tan /3 scenario in the pseudoscalar Higgs itla versus 
light chargino plane. In fig. [TS] the resulting change in the \ 2 is given for fixed tan/3 = 35. 



For small chargino masses there is a sharp increase in the \ 2 due to the corresponding mass 
limit, see above. The highest values for Rb can be obtained for small values of ijia and m;±, 



see fig. 16. As in the low tan/3 case the enhancement of Rb must be compensated by a 



decrease of a s , see fig. [17], and the change of R c is small, less than 0.0006 within the given 
parameter plane. The mixing angle, shown in fig. 18, is mainly determined by the Rb~, s -y 
rate, which can be fitted in the whole m^-m^i plane, see fig. |I9| . 

4.4 Best Solutions 

Standard Model Fits: 

The predictions of the SM are completely determined by the set of six input parameters 
M z , m t , the SM Higgs mass M h , a s , a em and G^. The error of the muon decay constant is 
so small that can be treated as a fixed parameter. The fine structure constant was taken 



to be l/a(M z ) = 128.89 ± 0.09||. The error on l/a(M z ) turns out not to be negligible: 



fixing l/a(M z ) underestimates the error on the Higgs mass by ~ 30%. The SM predictions 



were obtained from the ZFITTER package [12] and all the error correlations were taken from 



The fits were made with M z , m t , mjj, a s and a as free parameters, which resulted in 

M z = 91.186 ±0.002 GeV 

mt = 172.0^ GeV 

m H = 141+^f GeV 

a s (M z ) = 0.1197 ±0.0031 

l/a(M z ) = 128.90l8:™g 
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The minor differences to the results given in Q originate from the inclusion of the Rb-> s -y 
ratio. The electroweak mixing angle can be determined from these parameters: 



sin 2 Ojys = 0.2316 ± 0.0004. 



Both the strong coupling constant and the electroweak mixing angle sin 2 9j^g have been 
determined at the scale Mz in the MS renormalization scheme. With m t = 175 GeV 
sin 2 ^jg^- = sin 2 9 l ^Jl to an accuracy better than 0.0001. The quoted errors have been de- 
termined using MINOS fl~4j. Further details of the procedure are described elsewhere, see 



for example [37, 38]. The x 2 /d.o.f of the SM fit is 19.6/15 which corresponds to a prob- 
ability of 19%. Here, the main contributions to the x 2 originate from sin e l e f f from SLD 
(Ax 2 = 4.9), R b (Ax 2 = 3.1) and A b FB (A X 2 = 3.5). The prediction of R c is good. If R b ^ sy 
is excluded from the fit, one obtaines a x 2 /d.o.f of 18.9/14, corresponding to a probability 
of 17%, in agreement with ||. The correlation parameter between win and m t for the best 
fit is approximately 0.7; this strong correlation is shown in fig. ^t]. One observes that the 
upper limit on the Higgs mass is obtained for m t above 175 GeV; however, the upper limit 
is sensitive to sin 2 Q^fl as shown by the dashed contour in fig. |2D|, where the precise value 
of sin 2 from SLD was excluded from the fit. The dependence of sin 2 O^fl on the SM 
Higgs mass is approximately logarithmic (see fig. [pp. The LEP data without SLD yields 
ran = 24ll 2 23 GeV, while the SLD value of sin 2 6^ corresponds to m H = 16 ig 6 GeV; both 
rrih values are indicated in figure ETI (by the square and the circle, respectively). The SLD 



value is excluded by the lower limit of 58.4 GeV from the combined LEP experiments [[IT 
so more data is eagerly awaited. 

The Ax 2 dependence of the Higgs mass is shown in fig. ^2] for various conditions. 
Clearly, the sin 2 from SLD gives a large weight, while the new value from R b plays only 
in minor role in contrast to the previous value 0. 

MSSM Fits and Comparison with the SM: 

In order to obtain the best MSSM fits the assumption M 2 = 3|/i| is dropped and M 2 is 
treated as a free parameter. As discussed in section PTTJ the dominant contributions vary for 



the high and the low tan (5 scenario. The preferred tan (3 values for these scenarios are around 
35 and 1.6, respectively. Since the fit is not very sensitive to the precise tan/5 value, it was 
fixed to these values. The fitted MSSM parameters and the corresponding SUSY masses are 
given in table |3|; the predicted values of all observables and their pulls are summarized in 
table |2|. For the MSSM Mz, and a were treated as fixed parameters because of their small 
errors compared with the uncertainties from the other parameters. The MSSM prediction of 
the W-boson mass is always higher than the SM one, but the values of the strong coupling 
constant, the electroweak mixing angle and the top mass are very similar in the MSSM: 

a s (M z ) = 0.116 ±0.005 
m t = 172 ± 5 GeV 
sin 2 0^ = 0.2315 ±0.0004. 
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These values are for the low tan/3 scenario, but for high tan/5 the same values are obtained, 
except for a s = 0.119 ±0.005 in that case. The remaining parameters are given in table || A 
direct comparison to the SM fits is given in figs. [24|425|. The SM fit to the 20 measurements 
of table [| with 5 parameters yields x 2 /d.o.f . = 19.6/15 which corresponds to a probability 
of 19%, while the MSSM fits correspond to probabilities of 17% (tan/3 = 1.6, x 2 /d.o.f. = 
16.6/12) and 11% (tan/3=35, x 2 1 'd.o.f '. = 18.1/12). In counting the d.o.f the insensitive 
(and fixed) parameters are ignored, so only the parameters given in table |3| are considered. 
The difference in x 2 between SM and MSSM fits is mainly caused by Rb, which is better 
described in the MSSM, although the difference in x 2 is insufficient to distinguish between 
the models. 

The Higgs mass is not an independent parameter in the MSSM, since the couplings in the 
Higgs potential are gauge couplings, which limit the mass of the lightest Higgs to a rather 
narrow range |]39[|. The high tan/? needs a light pseudoscalar Higgs mass. As the lightest 
Higgs mass is strongly correlated with the pseudoscalar Higgs mass, it is also low. Similar 
Higgs values in the MSSM model were obtained in ref. |4(^ . 

An interesting point is the fitted value of a s (Mz). In previous analysis with the high 
values of Rb, a s (Mz) in the MSSM (« 0.11) was always significantly smaller than the 
SM value(0.123 ± 0.005[jllJ), which supported the low energy values from deep inelastic 
scattering (DIS) (0.112 ± 0.005|TTJ) and lattice calculations of the heavy quark splittings 
(0.110±0.006[[TTJ). However, the discrepancies between the low energy a s values and the LEP 
data have practically disappeared at the Warsaw Conference El[: the Standard Model value 



(0.120 ± 0.003, see above) is now in agreement with DIS measurements (0.115 ± 0.005 |f41~l), 
lattice calculations 0.117 ± 0.0030 and the world average 0.118 ± 0.0030. The MSSM 
values of a s are in good agreement with these other determinations, as shown in fig. 26. 

The particle spectrum for the best fits, as shown in table [5], suggests that some SUSY 
particles could be within reach of LEP II. Unfortunately, if the stop-, chargino- and/or Higgs 
mass are well above the discovery reach of LEP II, the x 2 °f the fit increases at most up to 
the SM value, since these particles basically decouple as soon as they become heavier than 
the heaviest SM mass, say the top mass of about 200 GeV. So one cannot get upper limits 
on these particles, since the probability changes only a few percent between the SM and 
MSSM. 



5 CMSSM and R b 

In the MSSM fits discussed above the lightest stop is mainly right-handed, while the left- 
handed stop has to be heavy. If both would be light, then all other squarks would likely be 
light, which would upset the good agreement between the SM and all other electroweak data. 
A large mass splitting in the stop sector needs a very artificial fine tuning of the few free 
parameters in the Constrained MSSM, which assumes unification of gauge and b-r Yukawa 
couplings ||38||. This is obvious from the sfermion mixing matrix for the stop quarks, eq. E3~ 



The D-terms proportional to cos 2(3 are negligible for tan/? fa 1. If one of the diagonal 
elements is much larger than m t , the off-diagonal terms of the order m t will not cause a 
mixing and the difference between the left- and right-handed stops has to come from the 
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evolution of the diagonal terms (for the notation see ref . [38] ) : 



dM% 16 I 

_9l = (_a 3 M 3 2 + 3a 2 M 2 2 + -aiM 2 ) (32) 

-[y t (M| 3 +M? 3 +m 2 ff2 +4 
r2 , 7i>r2 



+n(M4 + + + Aimi)} 

^ = (f^ + ^M?) ,33) 
-2y t (M| 3 + M? 3 +m^+^m 2 ) 

One observes that the difference between left- and right handed stops, denoted by Q3 and U3, 
respectively), depends on the Yukawa couplings for top and bottom (Y t , Y b ) and the trilinear 
couplings A t Q,y For low tan/5 Y b is negligible, while A t and Y t are not free parameters, 



since they go to fixed point solutions |$8fl. Therefore there is little freedom to adjust these 
parameters within the CMSSM in order to get a large splitting between the left- and right- 
handed stops. 

In addition, problems arise with electroweak symmetry breaking, since this requires the 



Higgs mixing parameter /i to be much heavier than the gaugino masses while R b re- 
quires low values of \i for a significant enhancement (since the chargino has to be preferably 
Higgsino-like). In conclusion, within the CMSSM an enhancement of R b above the SM is 
practically excluded. 



6 Conclusions 

Both the MSSM and SM provide a good description of all electroweak data. The best 
X 2 /d.o.f in the MSSM (SM) is 16.6/12 (19.6/15), which corresponds to a probability of 17% 
(19%). The lower \ 2 of the MSSM is mainly due to the better description of R b , but the fit 
requires an unnatural large splitting in the stop sector, as discussed in the previous section. 
Since the final analysis of most of the available LEP data is still in progress, one has to wait 
and see if the present preliminary value of R b will indeed stay above the SM prediction. 
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Symbol 


measurement 


best fit observables 






SM 


MSSM 


tan (3 and pull 




pull 


1.6 


pull 


35 


pull 


M 7 [GeV] 


91.1863 ± 0.0020 


91.1861 


0.08 


yl.lo63 




yl.lo63 




r z [GeVl 

^ L J 


2.4946 ± 0.0027 


2.4958 


-0.45 


O Ad AO. 


-U.UU 


z.4y4U 


n no 
U.ZZ 


a h [nb] 


41.508 ± 0.056 


41.468 


0.72 


A 1 A fi1 
41.401 


n Q/i 


A 1 A AO 
41.44y 


1 n^ 




20.778 ± 0.029 


20.755 


0.80 


zu. i oy 


n 19 


9n 779 
ZU. ( ( Z 


n 99 
o.zz 




0.0174 ± 0.0010 


n m fin 
U.U10U 


1.41 


n m R9 

O.OIUZ 


i 9n 


n ni R9 

O.OIUZ 


1 99 
l.ZZ 


Rb 


0.2178 ± 0.0011 


0.2158 


1.75 


0.2174 


0.38 


0.2168 


0.92 


R c 


0.1715 ± 0.0056 


0.1722 


-0.13 


0.1707 


0.14 


0.1708 


0.12 


A b 

FB 


0.0979 ± 0.0023 


0.1022 


-1.87 


0.1031 


-2.26 


0.1031 


-2.24 


A c 
FB 


0.0735 ± 0.0048 


0.0731 


0.10 


0.0736 


-0.02 


0.0736 


0.01 


A b 


0.863 ± 0.049 


0.933 


-1.45 


0.9353 


-1.49 


0.9356 


-1.50 




0.625 ± 0.084 


0.667 


-0.50 


0.6678 


-0.51 


0.668 


-0.51 




0.1401 ± 0.0067 


0.1460 


-0.88 


0.1470 


-1.03 


0.1466 


-0.97 




0.1382 ± 0.0076 


0.1460 


-1.03 


0.1470 


-1.16 


0.1469 


-1.14 


sin 2 9 l e e f pt f (Q FB ) 


0.2320 ± 0.0010 


0.2316 


0.35 


0.2315 


0.48 


0.2315 


0.46 


M w [GeV] 
1 - M^/Mf 


80.356 ± 0.125 


80.355 


0.01 


80.403 


-0.38 


80.428 


-0.58 


0.2244 ± 0.0042 


0.2235 


0.23 


0.2225 


0.45 


0.2220 


0.56 


m t [GeV] 
sm 2 9 l e f f (A LR ) 


175 ± 6. 


172.0 


0.50 


172.5 


0.42 


172.0 


0.50 


0.23061± 0.00047 


0.2316 


-2.21 


0.2315 


-1.94 


0.2315 


-1.97 


-Rfe^S7/10~ 4 


2.32 ± 0.67 ±0.5 


3.1 


-0.86 


2.43 


-0.12 


2.33 


0.0 


l/a{M z ) 


128.896 ± 0.09 


128.905 


-0.10 


128.89 




128.89 





Table 2: Measurements 0] and the predicted results of the fits with minimum % 2 . 
The pulls are defined by (measurement - predicted value) / error of the measurement. The 
second error for Rb^sj has been added to take care of the uncertainty by the renormaliza- 
tion scale used for the calculation of that quantity. For the MSSM fits Mz and \ ja[Mz) 
were taken as fixed parameters, because their uncertainties are negligible compared to un- 
certainties arising from the soft breaking parameters. Leaving them free does not change the 
results. 
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Fitted SUSY parameters and masses 


Symbol 


tan ,3=1.6 


tan/?=35 


m f [GeV] 


172 


172 


a s 


0.116 


0.1190 


M 2 [GeV] 


113 




//[GeV] 


60 


111 


m h [GeV] 


48 


187 




-0.18 


0.04 


m A [GeV} 




50 


Particle Spectrum 


m- tl [GeV] 


« 1 TeV 


m h [GeV] 


48 


187 


m g -[GeV] 


1 TeV 


m ; -[GeV] 


0.5 TeV 


m-±[GeV] 


149 


1504 


m^[GeV] 


84 


111 


m^o [GeV] 


54 


107 


m^o [GeV] 


64 


114 


m^o [GeV] 


100 


722 


m^o [GeV] 


150 


1504 


m h [GeV] 


109 


50 


m H [GeV] 


« 1.5 TeV 


112 


m A [GeV] 


1.5 TeV 


50 


m^±[GeV] 


« 1.5 TeV 


123 


X 2 /d.o.f. 


16.6/12 


18.1/12 


Probability 


17% 


11% 



Table 3: Values of the fitted parameters (upper part) and corresponding mass spectrum 
(lower part). On the right hand side the results of the optimization for high tan/5 are given. 
The dashes indicate irrelevant parameters which were chosen high. 
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Figure 4: Dependence of the chargino masses (solid line=lightest one) on the parameter \i 
for M 2 = \fi\ and M 2 = 3\/jl\ for tan/3 = 1.6, a s ~ 0.117 and m 4 - 2 « 60 GeV. The shaded 
regions indicate chargino masses less than 65 GeV which are excluded by LEP 1.5 and 
chargino masses less than 96 GeV, which is the discovery reach for LEP II. Note that two 
light charginos are easier obtained, if /J ~ M 2 . 




|0. (GeV) 



Figure 5: Dependence of the lightest chargino mass on the parameter /z for tan/3 = 35 
and M 2 =1500 GeV. The shaded regions as in fig. [|. The heavy chargino mass is close to 
M 2 = 1500 GeV(not indicated). 
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tan(3=1.6 




m-^GeV) (|^<-40 GeV) 



a 





\% 


\ \ \ \ 

>\ AA A- A 












1 i i i i i , 


50 52.5 


55 57.5 


60 62.5 65 67.5 70 72.5 

m-^GeV) (u>-40 GeV) 



Figure 6: Rb in the light stop versus light chargino plane with M 2 = 3|/x| and tan/3 = 1.6. 
The upper part shows the solution with < —40 GeV, in the lower part the one with 
ji > —40 GeV is displayed. In the latter solution quite high values for R b are possible, as 
can be seen in the figure. The dashed line in the upper plot indicates the old 2a lower limit 
of Rb- Recent updates of electroweak data yield Rb = 0.2178 ± 0.0011. 



tan[3=50 




50 55 60 65 70 

m-^GeV) 



Figure 7: Rb in the versus light chargino plane with M 2 = 3\/i\ for the high tan/3 
solution, fi was chosen positive here. In this case choosing the opposite sign for fi doesn't 
change Rb- 1fl 




Figure 8: Dependence of the absolute \ 2 f° r A 4 > — 40 GeV (left side) and // < —40 GeV 
(right side), using M 2 = 3|/j|. No optimization of parameters was performed, but they were 
fixed to values near the minimum. 




Figure 9: The same as |[ but for M 2 = \fi 
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Figure 10: The A% 2 in the light stop and light chargino plane for tan/3 = 1.6. At each 
point of the grid an optimization of m t , a s and the stop mixing angle 4> m i X was performed 
with /I > -40 GeV and M 2 = 3|/x|. 



21 




Figure 11: Rb in the light stop and light chargino plane. Optimization as in fig. [H| 
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Figure 12: a s in the light stop and light chargino plane. Optimization as in fig. |l(]. 
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-0.175 
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Figure 13: Stop mixing angle 0, 
determined by the Rb^ s ^ rate. Optimization as in fig. [10 



in the light stop and light chargino plane. It is mainly 
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Figure 14: i2&_» S7 in the light stop and light chargino plane for tan/3 = 1.6. For chargino 
masses higher than 60 GeV (and /i > 0) the predicted value is close to the CLEO measure- 
ment of 2.32 ± 0.67 x 10" 4 . Optimization as in fig. [K| 
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Figure 15: The A% 2 in the pseudoscalar Higgs and light chargino plane for tan/3 = 35. For 
each given tjia and light chargino mass an optimization of m t , a s , m t ~ 2 and the stop mixing 
angle 4> m ix was performed. The irrelevant parameter M 2 was set to 1500 GeV. 
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29 




Figure 19: -R&^ S7 in the pseudoscalar Higgs and light chargino plane for tan/3 = 35. The 
prediction is close to the CLEO measurement of 2.32±0.67x within the whole parameter 
space. Optimization as in fig. [L5. 
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Figure 20: A% 2 = 1 and A% 2 = 4 contour lines for all electroweak data including sin 2 O^Fj 
from SLD (continous line) and without it (dashed line). The stars indicate the best fits. 
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Figure 21: Dependence of the SM sin 2 d^Fj on the Higgs mass. The top mass m t = 
175 ± 9 GeV was varied within its error, as shown by the dashed band labelled SM (upper 
(lower) boundary m 4 =166(184) GeV). The SLD and the LEP measurement of sin 2 Q^f) are 
also shown as horizontal bands. The star and the square are respectivly the results of 
the combined fit to SLD and LEP data and LEP data without the SLD measurement of 
sin 2 Oeff, whereas the circle indicates the Higgs mass corresponding to the SLD measurement 
of sin 2 Clearly, the SLD value yields a Higgs mass less than the combined LEP limit of 
58.4 GeV (shaded area). 



io io 2 io 3 



m H (GeV) 



Figure 22: Dependence of the SM A% 2 on the Higgs mass for a free top mass, taking all 
data (continous line), all data without the SLD measurement of sin 2 O^Fj (dashed line) and 
all data without Rb (dotted line). 
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Figure 23: The A% 2 in the region of the best fit in the light stop and light chargino plane 
for tan/3 =1.6. Here the constraint on M 2 was dropped. At each point of the grid an 
optimization of m t , M 2 , a s and the stop mixing angle (j) m ix was performed with /i > —40. 
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Data / MSSM(tanB = 1.6) 

Data / SM 



0.85 0.9 0.95 



1 1.05 1.1 1.15 
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X /d.o.f = 19.6/15 
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0.85 0.9 0.95 



1 1.05 1.1 1.15 



Figure 24: Resulting observables for the fit given in table || for tan/? = 1.6. m~ b was fixed 
to 1000 GeV, rriA and the gluino mass were fixed to 1500 GeV. The remaining deviation of 
Rb from the SM prediction can be fully explained within the MSSM. 
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Data / SM 
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Figure 25: Resulting observables for the fit given in table |3] for tan/? = 35. m~ h was fixed 
to 1000 GeV, M 2 and the gluino mass were fixed to 1500 GeV. It is possible to improve the 
prediction of Rb with Supersymmetry even for high values of tan (3, but the result is not as 
good as for low values. 



GLS sum rule 
Bjorken sum rule 
QQ Lattice 
bb Lattice (prel.) 
deep inelastic 
Y decay 
X decay 

e + e" eventshape LEP 
SM 

MSSM tanf3=1.6 
MSSM tan(3=34 



-O 



-O 



-o- 



World Average 



o 
o 



o 



aJM 7 ) 



Figure 26: Comparison of different measurements of a s with our fit results, labeled SM and 



MSSM. The data has been taken from f 11 and 41 
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